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Abstract 

Properties of the Einstein equations in a coordinate gauge based on expand- 
ing null hypersurfaces foliated by metric 2-spheres are described. This null 
quasi-spherical (NQS) gauge leads to particularly simple analyses of the char- 
acteristic structure of the equations and of the propagation of gravitational 
shocks, and clarifies the geometry of timelike boundary condition. A feature 
of the NQS gauge is the use of the standard S ("edth") operator on S 2 to 
express angular derivatives, and the consequent use of spin- weighted spherical 

harmonic decompositions of the metric fields. 
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The objectives of numerical relativity in particular, force us to address the problem 
of constructing explicit formulations of the full Einstein equations. There are two main 
avenues of approach, respectively using either a 3+1 ADM representation []]], which leads to 
a Cauchy initial value problem, or a characteristic coordinate, an approach first considered 
byBondi g§. 

Both the second order and the recently popular first order [§-0] formulations of the 3+1 
equations involve coordinate gauge freedoms which must be eliminated through gauge-fixing 
conditions which typically involve subsidiary differential equations. First order formulations 
in particular introduce large numbers of auxiliary variables to parameterise the field, and 
these variables must satisfy additional compatibility relations. Whilst these features do 
not present difficulties of a conceptual or theoretical nature, they do cause difficulties in 
the practical numerical formulation of the equations — in disentangling geometric from 
coordinate gauge and numerical algorithm effects 0], in formulating the numerical algorithm 
to ensure that the constraints and compatibility relations are preserved by the evolution , 
and in simply constructing codes handling large number of variables. 

Characteristic methods ||, [TI| on the other hand typically involve very little coordinate 



freedom and lead to formulations of the equations involving a small number of variables, 
with few compatibility conditions. However, the advantage of simpler representations of the 
equations is balanced by the disadvantage that characteristic coordinates are generically not 
globally applicable, for either topological or caustic reasons. However, a characteristic-based 
coordinate system can be very useful when available — for example, in describing asymptotic 



structure p[|, or for describing discontinuity properties of vacuum spacetimes [1 1 



The Newman-Penrose formulation in terms of null tetrad based connection ("spin' : 



coefficients [|HJ can be used to describe the characteristic structure of the Einstein equations, 
either by regarding the equations as determining the propagation of the spin coefficients 
I3| , [l4|| , or by introducing coordinates and determining the spin coefficients in terms of 



the metric parameters (eg. |TtJ). The former ("pure NP") technique introduces numerous 
compatibility differential relationships amongst the spin coefficients, which suggests that 



more variables have been introduced than are necessary to fully describe the connection. 
The latter approach bridges the gap between coordinate- and connection-based techniques, 
and is the starting point for this paper. 

Ideally, a "good" parameterisation of the Einstein equations will have limited or no 
gauge freedom (to suppress unphysical gauge modes in the numerical evolution, for example 
g), with parameters free of constraints and having a direct relation to known radiation 
parameters, and leading to simple expressions for the Einstein equations, to facilitate analysis 
of the local existence, decay etc. 

The null quasi-spherical gauge presented here is a modification of the usual characteristic 
formulations ||, |K| which satisfies many of these criteria. It combines a characteristic 
gauge (exploiting the ode structure of the Einstein on a null hypersurface |J) with the 
reduced gauge- dependence of the null tetrad methods |12[], but avoids the complications 



of the spin coefficient integrability conditions. The quasi-spherical condition means that 
all derivatives can be presented using radial and time partial derivatives, combined with 
canonical differential operators on S 2 with the standard metric. 

An interesting feature is the central role played by the NQS metric field (3, which deter- 
mines in particular both the shear a^p of the outgoing null curves of the null coordinate, 
and the intrinsic geometry of the null hypersurfaces. Given (3 on one null hypersurface, 
the remaining NQS parameters are determined by integrating a system of radial ordinary 
differential equations which generalise the usual null Raychaudhuri equation (which is just 
the first equation (|j) of the system). Heuristically, (3 may be regarded as describing the 
ingoing gravitational radiation passing through the outgoing null hypersurfaces. 

The radial system of equations for the metric parameters is of course well-known in similar 
situations and can be understood in terms of the NP formulation of the equations for 

Einstein tensor terms $00, $01 > $02 and $n + 3A, which can be expressed solely in terms of 
derivatives tangent to the null hypersurfaces of the NP spin coefficients. 

However, closer examination of the explicit NQS coordinate representations shows that 
the hypersurface equations may be expressed in terms of hyper surf ace-tangential derivatives 



of the metric functions and just a subset of the connection parameters (see Eqns. (|8|-[lTD 
below). In all, we obtain a complete parameterisation of both the metric and the connection 
in terms of just 12 real variables, made up of 2 real and 2 complex variables for each of the 
metric (u,v,(3,j) and the connection (H, J,Q, K). Not only is this cheaper than the NP 
parameterisation (which requires 12 complex variables for the connection alone), but it is 
also explicit in the sense that there are no auxiliary compatibility conditions. The resulting 
expressions for the Einstein tensor (§-[ll]), (|i6Hi7|), (|T^) and Weyl tensor (p2Vp6|) components 
are comparable in simplicity with the usual NP expressions, yet have the advantage that 
they are completely explicit and involve fewer variables. 

The reduced set (H, J, K, Q) of connection parameters was found by careful examination 
of the explicit expressions for the hypersurface Einstein tensor components in the NQS 
metric. It seems quite likely that a similar examination using instead other characteristic- 
based metric forms (eg. Bondi, Newman- Unti) will uncover other reduced ("beyond NP") 
parameterisations for the connection. 

Because of the very explicit form of the Einstein equations in the NQS gauge, it enables 
us to give simple descriptions of several relatively well-known computations, such as: the 
null hypersurface matching conditions and the propagation laws for gravitational (vacuum) 
shock waves [flllJlql ; a description of the free boundary data for the Einstein equations on 
a timelike boundary surface; and the relationship between the outgoing Einstein tensor 
component and the geodesic parameterisation of the outgoing null rays. Two applications 
which are not discussed here, but which are quite feasible within the NQS framework, 



are the rederivation of the linearised Einstein equations |]16|-[T8|], and the asymptotic form 
of a vacuum NQS metric. The application which motivated the present study, namely the 
numerical evolution of of a vacuum black hole spacetime in the NQS gauge, will be described 



elsewhere |19|j, 
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I. NQS EINSTEIN EQUATIONS 



The central idea is to assume that the radial function foliates the null hypersurfaces 
by metric spheres, of area 4irr 2 . A similar foliation was exploited in to construct 
spatial initial data satisfying the Hamiltonian constraint. Denoting the null coordinate by z 
(cf. ||11|| ) and letting (p denote standard polar coordinates on the spheres, the most general 
compatible metric (assuming that r is a valid coordinate) is [] 

ds 2 N Q S = —2m dz{dr + v dz) + 2 \rQ + (3 dr + 7 dz\ 2 , (1) 

where we use a complex notation for simplicity, with 

9 = -±(d$ + i sin $ dtp), 

(3=^-i(3% 7 = ^(V-i 7 2 ). 

Thus u, v and 7 are metric parameters, with u, v real and u > 0, and we may regard 
/3, 7 either as spin-1 (complex) fields or as vector fields tangent to the spheres, via (3 ~ 
^dv + fi 2 esc i? d^. 

The metric form Eqn. ([]]) differs from other parameterisations of the metric with a 
characteristic coordinate 0, 0, |T2|, PH] , P^J principally in the choice of radial coordinate. 
For example, Bondi-Sachs use a luminosity parameter determined by the volume form in 
the angular directions [^3|] (with the consequent disadvantage that the luminosity coordinate 
depends strongly on the specific choice of angular coordinate labelling of the outgoing null 
curves) , and Newman-Unti use the geodesic affine parameter (and thus the radial coordinate 
is not determined uniquely from the intrinsic induced metric on the null hypersurfaces). 



Our metric has signature +2, the expressions below for the NP spin coefficients have been 
adjusted to give the usual signs, and the curvature tensor sign convention has been chosen so that 
Gu = 2 ( I ) oo is non-negative for physically reasonable matter. 
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However, the NQS radius and radial curves are determined by the metric S 2 foliation and 
standard polar coordinates, so that the NQS radial curves (z, (p) = const, do not in general 
coincide with the null generating curves. 

For the purposes of this paper it is sufficient to regard the metric form as an ansatz 
(assumption) rather than as a coordinate condition. The important question of whether this 
ansatz is generic, in the sense that coordinates satisfying the NQS conditions may be found 
in all metrics in some open subset of the space of "all" metrics, does not yet have a complete 
answer. Some of the evidence which suggests that the NQS coordinates are generically 
available is discussed below; we note also that the gauge used by Regge, Wheeler and Zerilli 
T6| , |T7|1 in the study of the linearisation of the Einstein equations about Schwarzschild, is just 



a linearised NQS gauge. Thus the NQS gauge could be used to extend the Regge- Wheeler- 
Zerilli perturbation analysis to include nonlinear effects. (The RWZ equations in the NQS 
gauge were rederived in fll8| .) 

The complex notation suggests encoding the angular derivatives using the differential 
operator <3 ("edth" [21|, defined on a spin-s field by 



^ \ dfi sin $ dip / 
All the standard differential operators on S 2 may be expressed in terms of 5, for example: 

div/5 = 3^ + 3/3 = /3.\ + /3. 2 2 
curl/? = i (8/3 - 30) = /3 1;2 - /3 2;1 
A0 = (58 + 88)0. 

Here the indices refer to the S^-orthonormal frame d$, esc $ d v , and the standard 5* 2 covariant 
derivative. The spin-weighted spherical harmonics Y s i m are then eigenfields of the Laplacian, 
satisfying AY s i m = (s 2 — 1(1 + l))Y s i m , \s\ < I, and Y s i m may be taken proportional to <S s Yi m 
for s > and d~ s Yi m for s < 0, where Y[ m are the usual spherical harmonic functions. 

To compute the curvature of ds 2 N Q S we introduce the vector frame (£, n, m, m), null with 
respect to dsj^gg. 
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£ = d r - r^d/s, 

n = u~ 1 {d z -r~ l d 1 -v(d r -r~ 1 df 3 )), (3) 

m = — -= (d$ — i esc tidy), 
ry 2 

and its associated dual (null) coframe dr + v dz, u dz, r 9 + (3 dr + 7 dz and r 9 + (3 dr + 7 dz. 
Note that denotes the angular directional derivative (3 l d$ + /3 2 csc^d v , which coincides 
with the S^-covariant derivative = /3<3 + /3S when acting on functions. 
It will be useful to introduce the differential operators 

V r := d r - r- l V p = d r - r~ l ((3Q + /35 s ), 
V z := d z - r-V 7 = d z - r- 1 ( 7 8 + 78), 

which have a natural interpretation as spin weight operators. Since we consider here 
only integer spins, we may use the complexified cotangent bundle T^ 1,0 ^ S 2 as a model for 
the spin-1 line bundle C 1 . From C 1 we construct the spin-s complex line bundle C s over 
R x M + x S 2 , whose sections are just spin-s fields depending also on the parameters (z,r). 
Then the action of TV, T> z on spin-s fields is defined via the standard covariant derivative 
on S 2 , extended naturally to T^ 1,0 'S 2 and C s . 

The structure of the hypersurface equations Q^-jTT[) motivates the introduction of the 
NQS connection parameters H, J, K, Q, Q ± , defined by: 

H = u _1 (2 -div/3), 
J = div7 + w(2-div/3), 

K = v(5f3 — 87, (4) 
Q = rV z f3 - rV rl + 7, 
Q ± = u - 1 {Q±du). 

These have well-defined spin-weights under S 2 frame rotations: u, v, H, J have spin (and 
are real), (3, 7, Q, Q ± have spin 1, and K has spin 2. 
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The connection V of ds 2 N q S may be described in terms of the NQS parameters either via 
the Cartan connection coefficients uj a b c = g(a, V c 6), a,b,c — £, n, m, m, or via the Newman- 
Penrose spin coefficients |[L2|| : 

KNP = ^>lml = 
&NP = ^Imm = r 

p NP = UJimrn = - \r~ 1 {2 - dlV/3) 

e^vp = |(^W + Urhmi) = \u~ lr D r u + i r~ x cot i? Im(/3) + |curl/?J 
\u^ x T) r v + ir _1 M _1 coti?Im(7 — + ^(curby — ucurl/3)j 



«iVP + PNP = TTVP = ^Inm = <^mn£ = | r 1 Q + 
&NP ~ PNP = ^mfhm = - COt 1? 

r^rp = ^ mn = ~r Q 

AjVP = ^mnm = T~ U~ K 
fi NP = UJfnnrn = — |r 1 M 1 J 

u NP = uj mnn = r^u^Bv 

As noted above, the NQS parameters (u, v, (3, 7) and (H, J, K, Q) form a more effi- 
cient and complete representation of the metric and connection than the NP coefficients 
ajyp, . . . ,t np , in part because they are defined in terms of a specific coordinate- adapted 
frame (£,n,m,fh). The "compatibility" conditions are also much simpler than those re- 
quired of the NP coefficients. Clearly, H, J, K, Q are determined uniquely from the metric 
parameters u, v, j3, 7. Conversely, if j3, H, J, K are given then the remaining metric parame- 
ters (u,v, r y) may be reconstructed as follows: u is found algebraically from div/3 and H; 7 
is found by solving 

C ^ '= 37 + 2^3 dlV7 = J 2^j3 ~ K (5) 
on each S 2 ; and finally v — ( J — div7)/(2 — div/3). Observe that the reconstruction process 
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is local to the null hypersurfaces, since T> z derivatives are not involved, and that Q is not 
used — it is used instead to determine df3/dz. 

If P = then Cp = Co = S, an elliptic operator on S 2 acting on spin-1 fields, which is 
surjective and has 6- dimensional kernel consisting of the I = 1 spin-1 spherical harmonics 
F llm , m = —1, 0, 1. Consequently if (3 is small (and we will show elsewhere that the pointwise 
bound 

< 3- 1/2 (2 - div/3) (6) 

is sufficient), then (H) is solvable for 7, uniquely if the I = 1 components of 7 are prescribed. 
The kernel of Cq may be identified with the infinitesimal Lorentz transformations and con- 
formal motions of S 2 — this suggests we may interpret this ambiguity in the reconstruction 
of 7 as a gauge degeneracy 

Observe that the hypersurfaces are expanding, p^p < 0, exactly when div/5 < 2, and 
this condition follows from (|6]). The deformation term 9/3/(2 — div/5) in Cp (|5|) is equal to 
— \<Jnp/ Pnp, which is an invariant of the null geometry (type I in the terminology of fll"5||). 

To understand the genericity of the NQS ansatz, consider the effect of an infinitesimal 
variation h a b of the null metric dsj^Qgl^ on a null hypersurface Af. A calculation shows that 
the metric S 2 condition is preserved (at the linearised level) if the accompanying coordinate 
deformation vector X = r((fh + (m) + rf£ satisfies 

where ho = h a i ) m a m b 1 h 2 = h a },m a m b , and 

_ divC +^h 
J ~ ~ 2 - div/3 ' 

Thus if the size condition (||) holds then the linearisation of the NQS gauge can be enforced 
by solving (0) for £; moreover, the existence of NQS coordinates in axially symmetric Bondi 
metrics was established in [fT8[] . Consequently it is plausible to regard the NQS gauge as a 
coordinate condition rather than an ansatz. Note also that the NQS parameters u, v, /3, 7 
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have the degrees of freedom (6 = 10 — 4) expected of a general metric after removing all 
coordinate degeneracies. 

The characteristic structure of the Einstein equations || becomes particularly clear when 
they are expressed in NQS coordinates using the connection quantities H, J,K,Q. Let 
G a b = R a b — \Rg a b be the Einstein tensor, normalised by Gu = £ a £ b G a b = 2$ 00 . Then Gu, 
Gi m , G£ n , G mm are given by the hyper surf ace equations H 

rD r # = (^drv/3 - 2 _ diy/3 j H, (8) 

rP r Q" = (9/3 - uH)Q- + Q~8/3 + 255/3 + - ifSw + 2r 2 G £m , (9) 

rX> r J = -(1 -div/3)J + w- \u\Q + \ 2 - |wdiv(Q + ) - ur 2 G £n , (10) 

rV r K = (§div/3 + i curl/3) K - |J3/3 + |«3Q + + ±w(Q + ) 2 + ±w 2 G mm . (11) 

These formulae provide the justification for introducing the specific forms (|]) for the NQS 
connection parameters (if, J, iT, Q). 

If j3 is known on a null hypersurface A4 then (|3|-|TT|) form a system of differential equations 
for fields on S 2 along the integral curves of the null generators i of Af z . Although the general 
ode structure is well known |J, the form of the RHS is considerably simpler than those found 
in other gauges f2|, [fLOU , p3| - Note that the NP equivalents of (PHTT|) are easily determined: 
for example, Eqn. @ is equivalent to the NP propagation identity (with k np = and 
Gu = 2$oo) 

DPnp = P 2 NP + Wnp\ 2 + 2Re(e NP )p NP + $ 00 , (12) 

which is just the null Raychaudhuri equation. 

The remaining Einstein equations may be handled via the second Bianchi identity |||9] . 

Lemma 1 Let M be an expanding null hypersurface with generating null vector i, and let 
T be a hypersurface which intersects each of the null curves of M exactly once. Suppose that 
F a b is a symmetric 2-tensor satisfying the conservation law F ab ,b = 0, and 
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Fee = Fe n = Fe m = F mm = 0, 
and such that F nn = F nm = on T . Then F a b = everywhere on Af ' . 

Proof: Writing out F ab ' h = in NP form, using Fee = Fe n = F lm = F mm = and the 



bicharacteristic condition n^p = 0, gives |12 



= Rep NP F mm , (13) 
F>e{F nm ) = (2p NP + Pnp — 2£Np)F nm + o~NpF nfh 

+ (tTjvp - T NP ) (14) 

De{F nn ) = 2Re(p NP — 2e NP )F nn + D m F nm 
- Rep NP F min + Re ((2 (3 NP + 2ti np 

^np ) F nm ) . (15) 



Since the expansion p^jp ^ by hypothesis, Eqn. (J13|) gives F m?fl = 0; eliminating F mm turns 
Eqn. flilf) into a linear homogeneous ode along the integral curves, which gives F nm = by 
the initial condition and uniqueness. Finally, Eqn. (|15"D now reduces to a linear homogeneous 
equation for F nn with only the zero solution satisfying the initial condition. QED. 

When applied to the Einstein tensor, F^ = G a b, the lemma shows that in order to 
show that the Einstein equations G a b = are everywhere satisfied, it suffices to find metric 
parameters satisfying only the hypersurface equations (|8|-|TT|) everywhere, and the boundary 
(subsidiary 0) equations G nn = G nm = on one hypersurface T transverse to the expanding 
null foliation. Observe that the signature of the hypersurface T is immaterial; all that is 
needed is that T be transverse to the outgoing foliation. 

The following explicit NQS expressions for G nm and G nn show that the boundary equa- 
tions may be used to constrain the initial data on T for the hypersurface equations ([SHlTl): 

rV z (J/u) = v 2 rV r (J/(uv)) + |(div7 - vdivfi) J/u 

+ 2u~ 1 \K\ 2 -\/ Q +v - Av + ur 2 G nn (16) 
r V Z Q + = (v rV r + J - vftfi + d^)Q + - KQ + + 

+ 2u~ l rV r (udv) - (2 + icurij3)9v + 5 J - 2u~ 1 Jftu - 2ur 2 G nm . (17) 

11 



The terms on the right hand sides of (fL6| , |P7l) are determined on a single null hypersurface from 
the hypersurface equations and a gauge choice for 7/=i, thereby determining the evolution 
terms d(J/u)/dz and d(Q + )/dz. This constrains (J/u) and Q + on T - clearly when T 
is taken to be a level set of r (so d z is tangent to T), whilst for more general transverse 
hypersurfaces T, using knowledge of d(J/u)/dr and d{Q + )/dr. 

The remaining Einstein tensor equation G mm = is called the trivial equation in || by 
virtue of (1131) . In the NQS parameterisation, G mm is given explicitly by 

ur 2 G mm = rV r J - ±div/3 J - u\Q + \ 2 + |wdivQ + + KQ(3 + KB(3 + Q+Su + Q + 8u 

+ r 2 (vV r - V z ){vT x V r u) + u -1 r 2 Z> r (uZ> r v). (18) 



II. SOME APPLICATIONS 

The simplicity of the explicit NQS expressions for the Einstein tensor suggests many 
possible applications, and in the following we outline three: the formal structure of the null- 
timelike boundary value problem B ; matching conditions across a null hypersurface and the 



propagation of impulsive waves jlljl ; and the relation between the null intrinsic geometry 
(type I |[l5|) and the class of affine null geodesic parameterisations. 

A formal solution algorithm for the Einstein equations in the NQS gauge starts with /3 
defined on a null hypersurface Af z . Specifying H, J, K, Q~ on the intersection Tr\J\f z with a 
transverse hypersurface gives boundary (initial) conditions for the hypersurface equations, 
with [3\j\f z providing source terms. Solving the hypersurface equations for H, J, K, Q~ on 
M z and imposing a gauge fixing condition on the I = 1 components of 7 (eg. ji=\ = 0) 
then determines the remaining metric parameters u, v, 7. The definition Eqn. (|j) of Q~ 
determines d(3 /dz on Af z , which is used to evolve (3 to the "next" null hypersurface Af z+ Az- 
The boundary equations (|TB],[T7|) determine new boundary data for J/u, Q + , which are 
combined with arbitrary boundary data for u, K on T and the new seed field (3\n z+Az to 
repeat the process. Practical experience with a direct numerical implementation of this 
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procedure has been very good, which suggests that it might be possible to prove existence 
theorems by this approach. 

The freedom in the choice of u on T H M z should be interpreted as the freedom in 
choosing the outgoing null hypersurfaces. This can be seen as follows: at the linearised 
level, the coordinate tangent vector d z points from one null hypersurface to the "next" - 
the ^-translations are diffeomorphisms preserving the null foliation, and generated by the 
vector field d z . Adding a component tangential to the null hypersurfaces does not affect this 
property. Since the NQS null vector I = d r — r~ 1 (3 depends only on the intrinsic (type I) 
geometry of J\f z , the equivalence class of spacetime evolution vectors W preserving a given 
null foliation may be characterised by the inner product g(W,£). But for the deformation 
vector W = d z , this inner product is determined by u since u = —g(d z , £). Hence the choice 
of u > on T determines the intersection 2-surfaces M z D T, which in turn generate the null 
hypersurfaces Af z ; in other words, the boundary freedom in u is a gauge (NQS coordinate) 
freedom. 

The boundary data for K (essentially the incoming shear Xnp) is unconstrained and 
may be considered as representing the outgoing gravitational radiation injected into the 
spacetime from the region beyond the boundary T. In a similar fashion, we may regard the 
shear vector (3 (which is arbitrarily specifiable on one outgoing surface Af z ) as measuring the 
incoming radiation. 

Since (3 acts as a source for the remaining fields on the null hypersurface, we might expect 
NQS metrics with j3 identically zero to play a distinguished role. Such metrics admit a null 
geodesic congruence which is expanding, twist free and shear free, and thus correspond to 
the Robinson- Trautman spacetimes ||26|| . The NQS parameterisation of Robinson- Trautman 
metrics with S 2 cross-sections is described in E7[ and the remaining NQS gauge freedoms 



are classified in 28 



NQS coordinates may also be used to describe the junction conditions for null hyper- 
surfaces. The NQS approach has the advantage of simplicity, when compared with general 
coordinate methods []ll~l , |2"9| which must separately consider the effects of coordinate changes 
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on the discontinuity fields. On the other hand, the reliance on a special coordinate system 
apparently limits the applicability of the NQS analysis. 

Let M be two spacetimes with NQS coordinates and having null boundary pieces Af^ 
(one past, the other future), together with an identification map Af~ ~ Af + — Af . We 
use the notation [■] to denote the discontinuity of parameters across the matching surface, 
eg. [P] = (3^+ — Pm-. Supposing that the matching is an isometry along Af, ie. [P] = (this 
is the most natural matching criterion), we consider the effect of requiring that the Einstein 
equations be satisfied in the weak (distributional) sense. 

Let T be a hypersurface transverse to Af in M — M + U M~ . Since we may regard 
u\r as gauge, it is reasonable to require also that [u] — on T fl Af. Likewise we assume 
that the / = 1 spherical harmonic components of 7 are gauged to zero. The boundary 
equations (|I6| ), (|P7D show that [J/u] = 0, [Q + ] = on TflAf since their right hand sides are 
bounded on Af, from which it follows (using the definitions (|)) that [H] = [J] = [Q~] = 
on T fl Af. Considering the jump in (g) across Af shows that [H] satisfies the ordinary 
differential equation 

rV r [H] = (^drv/3 - - _ — j [H], (19) 

along the null generating curves of Af, and uniqueness with the initial data [H] = on TtlAf 
shows that [H] = on Af. Since [div/3] = we have [u] = also. Similar arguments using 
Eqn. © show that [Q~] = (and thus [Q\ = [Q + ] = 0) on TV, and using Eqn. (0), that 
[J] = 0. 

However, the discontinuity [K] on TC\Af is not constrained in similar fashion by gauge or 
boundary equation considerations, and thus Eqn. (|TTD yields a formula for the propagation 
along Af: 

rV r [K] = (|div/5 + i curl/5) [K]. (20) 

Using the gauge condition [7^=1 = on Af to invert Cp, from Eqn. (|5|) we recover the 
discontinuities across Af of the remaining metric coefficients, 
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and we find [d(3/dz] 7^ 0, using [Q] = and Eqn. (|[). From the NQS formulas for the Weyl 
curvature spinors 

rV'fo = (rX> r + 1 - 23/3)(3/3/u), (22) 
4r 2 ^x = (rl? r - 9/3)Q- - 3div/3 

- (3Q + gu)u _1 3/3 + 285/5, (23) 
2r 2 ^ 2 = ±r 2 (2G fa + G mfh ) - (1 - ±ff J) 

- 2u- 1 Rdp - |(gg + - 8g+), (24) 

4rV^ 3 = w(rZ> r - 3/3)Q + - u(rV z - dj)Q+ 
+ 2(rV r - 1 - ±i curl/3) (udv) 

- udJ - (3Q - 6u)K - 2uBK, (25) 
r 2 ^ = t> 2 (r£> r + 1 - 2dj3)(K/(uv)) 

- (rV z - 257) (K/u) + Q + dv + 55%, (26) 

we see easily that \l/o, ^1, ^2 and ^3 will be bounded, but ^4 will have a 5-function compo- 
nent —r~ 1 u~ 2 [K]S(z), giving an impulsive gravitational wave propagating along TV. 

Although the resulting metric is not continuous in NQS coordinates, there is a (non-C 1 ) 
coordinate change which makes the metric continuous, since the boundary identification is 



an isometry f30fl . There are two interesting points about this NQS construction: it provides a 
geometric coordinate condition which does not detect the "optimal" metric regularity (C°) , 
and secondly, it provides strong evidence for the existence of vacuum spacetime metrics 
admitting particular coordinates in which the metric is not continuous, but such that the 
full curvature tensor, computed in those coordinates, is well defined as a distribution. Such 
discontinuities are well beyond the reach of the usual Sobolev space approaches to local 
existence for the Cauchy problem for the Einstein equations, for which g £ if 2 seems to be 
the best plausibly possible (and g G if 5 / 2+e is the best proven result at this time). 
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The difference between characteristic and 3 + 1 approaches is again underscored by 
noting that whereas the Einstein equations across a null matching hypersurface may be 
satisfied with unbounded curvature components, determined by the shock transport law (pop , 
imposing the Einstein equations across a smooth spacelike matching hypersurface forces the 
metric to be C 2 ~ in Gaussian coordinates P§|] . 

From the dual reformulations of Eqs. (|lOl), (P]), 

rV z {uH) = rV r (uvH) + ur 2 G ln 

- u(l + Hv- HJ) + \u{\Q-\ 2 - divQ-), 
rV z (dp) = rV r (vd(3) - \ur 2 G mm 

+ (—v + i curhy — i v curl/3) 9/3 

+ \{uHK + J9/3) - \u{Q~) 2 + \udQ-, 

which correspond to interchanging the role of the inward and outward null directions, we see 
that the conditions div[/3] = 0, B[f3] = are required by the Einstein equations. Thus the 
isometry assumption [/?] = amounts to a constraint on the remaining, purely rotational, 
component of [/?]. 

Because the NQS foliation is intrinsic to the null hypersurfaces (determined by the type 
I geometry, in the terminology of Hl5|), the interaction between the Einstein tensor and the 
null geodesic parameterisation may be simply explained. It is easy to see that u~ x l is a 
null geodesic vector field, and thus the geodesic parameterisation is determined by the null 
type I geometry (ie. (3) and u, which is in turn determined from Gu = 2$oo by Eqn. (H), up 
to a choice of initial condition which amounts to fixing the scaling of the affme parameter. 
Thus, the effect of matter (as described by Gu > 0) is to slow down the null geodesic 
parameterisation, when compared to the parameterisation on an isometric null hypersurface 
in vacuum. This leads to a redshift of outgoing photons, as seen by an observer at constant 
NQS radius. 
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